Quantum Arnol'd diffusion in a rippled waveguide 
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We study the quantum Arnol'd diffusion for a particle moving in a quasi- 

ID waveguide bounded by a periodically rippled surface, in the presence of the 
time-periodic electric field. It was found that in a deep semiclassical region 
the diffusion-like motion occurs for a particle in the region corresponding 
to a stochastic layer surrounding the coupling resonance. The rate of the 
quantum diffusion turns out to be less than the corresponding classical one, 
thus indicating the influence of quantum coherent effects. Another result is 
that even in the case when such a diffusion is possible, it terminates in time 
due to the mechanism similar to that of the dynamical localization. The 
quantum Arnol'd diffusion represents a new type of quantum dynamics, and 
may be experimentally observed in measurements of a conductivity of low- 
dimensional mesoscopic structures. 
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As is well known, one of the mechanisms of the dynamical chaos in the Hamiltonian 
systems is due to the interaction between nonlinear resonances [1]. When the interaction 
is strong, this leads to the so-called global chaos which is characterized by a chaotic region 
spanned over the whole phase space of a system, although large isolated islands of stabihty 
may persist. For a weak interaction, the chaotic motion occurs only in the vicinity of 
separatrices of the resonances, in accordance with the Kolmogorov-Arnol'd-Moser (KAM) 
theory (see, for example, Ref. [2]). In the case of two degrees of freedom {N — 2), the 
passage of a trajectory from one stochastic region to another is blocked by KAM surfaces. 

The situation changes drastically in many-dimensional {N > 2) systems for which the 
KAM surfaces no longer separate stochastic regions surrounding different resonances, and 
chaotic layers of the destroyed separatrices form a stochastic web that can cover the whole 
phase space. Thus, if the trajectory starts inside the stochastic web, it can diffuse through- 
out the phase space. This weak diffusion along stochastic webs was predicted by Arnol'd in 
1964 [3], and since that time it is known as a very peculiar phenomenon, however, univer- 
sal for many-dimensional nonlinear Hamiltonian systems (see, for example, review [4] and 
references therein). 

Recently, much attention has been paid to the chaotic dynamics of a particle in a rip- 
pled channel (see, for example, Refs. [5-9]). The main interest was in the quantum-classical 
correspondence for the conditions of a strong chaos. Specifically, in Ref. [5] the transport 
properties of the channel in a ballistic regime were under study. Energy band structure, 
the structure of eigenfunctions and density of states have been calculated in [6,7]. The 
quantum states in the channel with rough boundaries, as well as the phenomena of quan- 
tum localization have been analyzed in [9]. The influence of an external magnetic fleld for 
narrow channels was investigated in [8]. These studies may have a direct relevance to the 
experiments with periodically modulated conducting channels. In this connection one can 
mention the investigation [10] of transport properties of a mesoscopic structure (sequence 
of quantum dots) with the periodic potential formed by metallic gates. 

In contrast with the previous studies, below we address the regime of a weak quantum 
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chaos which occurs along the nonhnear resonances in the presence of an external periodic 
electric field. Our goal is to study the properties of the quantum Arnol'd diffusion which 
may be observed experimentally. The approach we use is based on the theory developed in 
Ref. [11] by making use of a simple model of two coupling nonlinear oscillators, one of which 
is driven by two-frequency external field. 

We study the Arnol'd diffusion in a periodic quasi-one dimensional waveguide with the 
upper profile given in dimensionless variables by the function y = d + a cos x. Here x and y 
are the longitudinal and transverse coordinates, d is the average width, and a is the ripple 
amplitude. The low profile is assumed to be fiat, y = 0. The nonlinear resonances arise due 
to the coupling between two degrees of freedom, with the following resonance conditions, 

rj = T^/Ty = Uy/U^ (1) 

Here Ty is the period of a transverse oscillation inside the channel, Tx is the time of fiight 
of a particle over one period of the waveguide, and ujy are the corresponding frequencies, 
and 7] is the rational number. 




FIG. 1. Some of the coupling resonances for different values of r/, with the isoenergetic curve 
E = 160000 in the frequency plane. The locations of resonances are shown by dashed lines. 
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The mechanism of the classical Arnol'd diffusion in this system is illustrated in Fig. 1. 
Here some of the resonance lines for different values of rj are shown on the cu^, a; j,- plane. 
Also, the curve of a constant energy is shown determined by the equation 

ul+(^^J = 2mE. (2) 

Here E and m are the dimensionless kinetic energy and particle mass, respectively, which 
we set to unity in what follows. The neighboring coupling resonances are isolated one from 
another by the KAM-surfaces, therefore, for a weak perturbation the transition between their 
stochastic layers is forbidden. Such a transition could occur in the case of the resonance 
overlap, i.e. in the case of the global chaos only. In the absence of an external field, the 
passage of a trajectory along any stochastic layer (this direction is shown at Fig. 1 by two 
arrows) is also impossible because of the energy conservation. However, the external time- 
periodic field removes the latter restriction, and a slow diffusion along stochastic layers 
becomes possible. 

An external electric field with the corresponding potential V{y,t) = —foy{cosflit + 
cosft2t) gives rise to two main resonances with the frequencies cUy — Jli and cUy — ^12- Their 
locations are shown in Fig. 1 by the dashed lines. In order to calculate the diffusion rate, we 
consider a part of the Arnol'd stochastic web created by three resonances, namely, by the 
coupling resonance and two guiding resonances with frequencies fli and fl2- Therefore, we 
chose the initial conditions inside the stochastic layer of the couphng resonance. To avoid 
the overlapping of the resonances, we assume the relation fo/a = 1000 is fulfilled. 

For the further analysis it is convenient to pass to the curvilinear coordinates x[ in which 
both boundaries are flat [12]. The covariant coordinate representation of the Schrodinger 
equation has the following form, 

-2^1:^^-^ = 

where gik is the metric tensor, g = det{gij). Here we use the units in which the Plank's 
constant and effective mass are equal to unity. As a result, the new coordinates are 



where e — a/d. In these coordinates the boundary conditions are ■0(x',O) = ip{x',d) — 
and the metric tensor is 

/ -j^ ex' sin x' \ 



1+ecosi;' (l+ecosa;')2 

with the orthonormahty condition, 

/ i^ti^i^dS' = 5,,-. (6) 

If the ripple amphtude a is small compared to the channel width d, than keeping only the 
first-order terms in e in the Schrodinger equation (3), we obtain the following Hamilto- 
nian [12], 

H^Ho{x,y) + U{x,y), (7) 

where 
and 

A e / ^ . d 1 . d\ 

U — - [ 2 cos a;— — r- — 2y sm x——- y cosx— cos x — sm x— . (9) 

2 y oy^ oxay oy 2 ox J 

Here and below we omitted primes in coordinates x' and y'. 

Since the Hamiltonian is periodic in the longitudinal coordinate x, the eigenstates are 
Bloch states. This allows us to write the solution of the Schrodinger equation in the form 
ip{x, y) = exp{ikx)^k{x, y) where $fc(a; + 27r, y) = y). For an infinite periodic channel 
the Bloch wave vector k has continuous values, in particulary, —1/2 < k < 1/2 in the first 
Brillouin zone. 

By expanding ^k{x,y) in the double Fourier series the eigenstates can be written as 

i^\x,y) = e^'^j:^L€mi^,y) (10) 

n,m 



where 

^L(^,y) = /5e-^sin(^) (11) 

are the eigenstates of the unperturbed Hamiltonian HQ{x,y) with the corresponding eigen- 
values 

^nm=l{in + kr+^^ (12) 

for the considered case with L » 27r. Now we proceed to solve the system of algebraic 

equations for the coefficients c^^, 



1 / 7rW\ 

/ n' ,m' 



Here the matrix elements are 



a 
'2d 



^ {^n',n+l + ^n',n-l) Sm,m' + ''"l^UT' ((1 + 2(A; + u)) + (1 - 2{k + u)) (5„,,„_i) 

(14) 



Following Refs. [11], we analyze the dynamics in the vicinity of the main coupling resonance 
77 = 1 which is determined by the condition where LOn^ — Eno+i{k) — En^ik) — 

k + riQ + 1/2 and ujmo = -^mo+i ~ -^mo = '7r^(2mo + l)/2d?. In a deep semiclassical region 
where no ^ 1 and mo ^ 1, one can write a;^,, ~ no and ujjno ~ vr^mo/o?^. It should be noted 
that the similar resonance condition can be satisfied for negative Uq also in the case when 
—no ~ 'n-'^mQ/(P, which corresponds to the particles moving in the opposite direction. Below 
we use the fact that when |no| ^ 1 the two resonances (two sets of states), corresponding to 
no > and no < are not coupled for all k apart from the center and edges of the Brillouin 
zone. This fact is due to the anti-unitary symmetry of the Hamiltonian for all values of k, 
apart from the indicated above. The properties of the energy spectra and eigenstates for 
specific values k — ±1/2 and k — will be discussed elsewhere. 

In the vicinity of the resonance it is convenient to introduce new indexes r = n — no and 
p — r + {m — mo). Then, instead of the system (13) we obtain, 
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E{k)c^p = pUmoCrp +^[r^ + -p(P-^)j^rp+Yl Uk+no-r,p-r+mo,k+no-r' ,p' -r'+moCr'p' ■ (15) 

\ / r',p' 

Here we count energy E{k) from the level E^^^^^^{k). 

Numerical calculation of the solution of these equation gave the following results. First, as 
in Ref. [11], the energy spectrum consists of a number of Mathieu-like groups corresponding 
to the coupling resonance. These groups are separated one from another by the energy Uno- 
The structure of energy spectrum in each group is typical for a quantum nonlinear resonance. 
Inside the resonance, the lowest levels are practically equidistant, the accumulation point 
corresponds to the classical separatrix and all states arc non- degenerate. The states above 
the separatrix are quasi-degenerate due to to the rotation in opposite directions. 

In accordance with the spectrum structure it is convenient to characterize the states at 
coupling resonance by two indexes: group number q and s — level number inside the group. 
Correspondingly, the energy of each group can be written as follows, 

E,,,{k)^u;no{k)q + E^^{k), (16) 

where E^^ is the Mathieu-like spectrum for one group. The indexes q and s correspond to 
fast and slow variables characterizing the motion inside the classical coupling resonance. 

Let us consider now the dynamics of a charged particle in the rippled channel in 
the presence of the time dependent electric field described by the potential V{y, t) ~ 
—foy{cosQit + cosQ2^)- We assume that the frequencies Qi and Q2 are chosen to ful- 
fill the condition a;^^ = {Qi + (^2) /2 in order to provide equal driving forces for a par- 
ticle inside the stochastic layer of the separatrix under consideration. Specifically, we 
take, LOno — 400, Qi = 350, Q2 = 450, therefore, the period T of the perturbation is, 
T = 7- 27v/Qi = 9 • 27v/n2 ~ 0.126. 

Since the total Hamiltonian is periodic in time, one can write the solution of the non- 
stationary Schrodinger equation as ip{x,y,t) — exp {—iegt) UQ{x,y,t), where UQ{x,y,t) is 
the quasienergy (QE) function and Sg is the quasienergy. As is known, the QE functions 
are the eigenfunctions of the evolution operator U (T) of the system for one period of the 
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perturbation. The procedure to determine this operator was described in details in Ref. [11]. 
The matrix elements Uq^s,q',s'{T) of the evolution operator can be calculated by means of the 
numerical solution of the non- stationary Schrodinger equation. Then, the evolution matrix 
Uq,s,q',s'{NT) for periods can be easily obtained. 

Our goal is to analyze the dynamics of a particle placed inside the separatrix under the 
condition that the coupling and two driving resonances do not overlap. The evolution of 
any initial state can be computed using the evolution matrix as follows, 

C,,s{NT) = J2 Uq,s,q',ANT)C,,,s'{0). (17) 
q',s' 




FIG. 2. Time dependence of the variance Ag for different initial states at coupling resonance for 
the group with q = 0: (a) the lowest level with s = 0, (b) the above-separatrix level with s = 45, 
(c) the near separatrix level with s = 22. Here a = 0.01 and /o = 10. 



In Fig. 2 typical dependencies of the variance y^Hj = uo^^^Aq of the energy are shown 
versus the time measured in the number of periods of the external perturbation, for 
different initial conditions. The quantity Ag is defined as follows, 

A, = - E KA\ g = E ic,,.r- (18) 



The data clearly demonstrate a different character of the evolution of the system in 
dependence on the initial state. For the state taken from the center of the coupling resonance, 
as well as above the separatrix, the variance oscillates in time, in contrast with the state 
taken from inside the separatrix. In the latter case, after a short time the variance of the 
energy increases linearly in time, thus manifesting a diffusion-like spread of the wave packet. 

In order to characterize the speed of the diffusion, we have calculated the classical diffu- 
sion coefficient for the coupling resonance 77 = 1, in the comparison with the quantum 
one, Dq^ for different goffer amplitude and initial wave number k 7^ 0;±l/2, see Fig. 3. 
It was found that the quantum Arnol'd diffusion roughly corresponds to the classical one. 
However, one can see a systematic deviation which indicates that the quantum diffusion is 
weaker than the classical Arnol'd diffusion. This deviation is due to the influence of co- 
herence quantum effects which can be very strong even in a deep semiclassical region. As 
was shown in Ref. [13], these stabilizing quantum effects are enhanced for the motion inside 
narrow stochastic layers surrounding the nonlinear resonances. This effect is due a relatively 
small number Mg of quantum states belonging to the region of the separatrix, the fact which 
is crucial in the study of the quantum-classical correspondence for the systems with a weak 
chaos in the classical limit. 

logioD 
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FIG. 3. Classical Dd versus quantum Dq diffusion coefficients in the dependence on the ampli- 
tude a of the rippled profile. 



In this connection it is instructive to estimate the number of the energy states that 
occupy the separatrix layer in our calculations. We have found that for a = 0.01 of the am- 
plitude characterizing the upper profile of the waveguide, the number of stationary states in 
the separatrix chaotic layer is about Mg ~ 10, therefore, one can treat our results semiclas- 
sically. Our additional calculation has shown that with the decrease of the amplitude a, the 
number Mg decreases, and for l/s/a ~ 20 this number is of the order one. For this reason 
the last right point in Fig. 3 corresponds to the situation when the chaotic motion along 
the coupling resonance is completely suppressed by quantum effects. This effect is known 
as the "Shuryak border" [13]) which establishes the conditions for a complete suppression 
of classical chaos. 

Since the diffusion along the coupling resonance is effectively one-dimensional, one can 
expect an Anderson-like localization which is known to occur in the presence of any weak 
disorder in low-dimensional structures. Indeed, the variance of the QE eigenstates of the 
evolution operator is finite in the g-space. This means that eigenstates are localized, and the 
wave packet dynamics in this direction has to reveal the saturation of the diffusion. More 
specifically, one can expect that the linear increase of the variance of the energy ceases after 
some characteristic time. 

In order to observe the dynamical localization in our model (along the coupling resonance 
inside the separatrix layer), one needs to analyze a long-time dynamics of wave packets. 
Our numerical study for large times (see curve (c) in Fig. 2) have revealed that after some 
time t ~ 200T, the diffusion-hke evolution terminates for all range of the amphtude a. 
Specifically, on a large time scale the variance starts to oscillate around a mean value. This 
effect is of the same origin as the so-called dynamical localization which was discovered in the 
kicked rotor model [14,15]. One should note that the dynamical localization is, in principle, 
different from the Anderson localization, since the latter occurs for the models with random 
potentials. In contrast, our model is the dynamical one (without any randomness), and a 
kind of pseudo-randomness is due to the quantum chaos mechanism. 

In conclusion, we have studied the quantum Arnol'd diffusion for a particle in the quasi- 
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ID rippled waveguide under the time-periodic electric field. This diffusion is known to 
occur in the corresponding classical systems, below the threshold of an overlap of nonlinear 
resonances resulting in a strong chaos. As is known, in this case the chaotic motion is 
possible only for the trajectories inside narrow stochastic layers surrounding the nonlinear 
resonances. The classical Arnol'd diffusion is exponentially weak, however, it leads to the 
unbounded motion inside the resonance web created by the resonances of different orders. 

The problem under study was to find out whether the Arnol'd diffusion is feasible in 
quantum systems describing the electron motion in quasi- ID waveguides. Our results clearly 
demonstrate the peculiarities of the quantum Arnol'd diffusion and establish the conditions 
under which it can be observed experimentally. One of the most important results is that 
the quantum diffusion is typically weaker than the corresponding classical one. The analysis 
has shown that the reason is the influence of quantum coherent effects which are essen- 
tially strong for the motion in the vicinity of the nonlinear resonances, in the case when 
these resonances do not overlap. Moreover, we found that with a decrease of the perturba- 
tion responsible for the creation of stochastic layers, the quantum Arnol'd diffusion can be 
completely suppressed, thus leading to the absence of any chaos in the system. 

Another important result is the observation of the dynamical localization of the Arnol'd 
diffusion, which is manifested by the termination of the diffusion on a large time scale. 
We have briefly discussed the mechanism of this phenomenon, relating it with the famous 
Anderson localization known to occur in low-dimensional disordered structures. Thus, in 
addition to the Shuryak border, the dynamical localization is a new mechanism destroying 
the quantum Arnol'd diffusion. 

Recently, the weak electron diffusion was observed experimentally [16] in superlattices 
with stationary electric and magnetic fields. At certain voltage, the unbounded electron 
motion along a stochastic web changes the conductivity of the system and results in a 
large increase of the current flow through the superlattice. Similarly, one can expect that 
electron Arnol'd diffusion inside resonance stochastic layers can increase the high frequency 
conductivity of a rippled channel, and can be detected experimentally. 
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